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Abstract 

If neutrinos are to play a relevant cosmological role, they must be essentially 
degenerate with a mass matrix of the bimaximal mixing type. We study this 
scenario in the MSSM framework, finding that if neutrino masses are produced 
by a see-saw mechanism, the radiative corrections give rise to mass splittings 
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MSW) solar neutrino oscillations. This provides a natural origin for the Am^^i <C 
Am^j^ hierarchy. On the other hand, the vacuum oscillation solution to the 
solar neutrino problem is always excluded. We discuss also in the SUSY scenario 
other possible effects of radiative corrections involving the new neutrino Yukawa 
couplings, including implications for triviality limits on the Majorana mass, the 
infrared fixed point value of the top Yukawa coupling, and gauge coupling and 
bottom-tau unification. 
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1 Introduction 



If neutrinos are to play a relevant cosmological role, their masses should be 0{eV). 
In that case, since atmospheric and solar neutrino anomalies ^ indicate that mass- 
squared splittings are at most 10~^ eV^, neutrinos must be almost degenerate On 
the other hand, supersymmetry (SUSY) is a key ingredient in most of the extensions 
of the Standard Model (SM) which are candidate for a more fundamental theory. 
In this paper we will analyze, within the supersymmetric framework, under which 
circumstances the "observed" mass splittings between quasi-degenerate neutrinos arise 
naturally (or not), as a radiative effect, in agreement with all the available experimental 
data. 

This problem has been also addressed in a recent paper by Ellis and Lola 0, in 
which they treat the neutrino mass matrix, Aii,, as an effective operator, emerging at 
some scale. A, with the bimaximal mixing form. Then the renormalization group (RG) 
analysis shows that the splittings and mixings at low energy are not in agreement with 
observations. Here we take a more general point of view. Besides exploring the effective 
operator scenario, we focus our attention in the (well motivated) case in which this 
operator is produced by a see-saw mechanism^. This introduces crucial differences in 
the analysis. In particular, the form of Ai^, is modified by a first stage of RG running of 
the neutrino Dirac- Yukawa matrix and the right-handed neutrino mass matrix from the 
high energy scale (say Mp or Mqut) to A, which is identified with the Majorana mass 
scale. As we will see, this modification allows in many cases to reconcile the scenario 
with experiment, providing also a natural origin for the "observed" solar-atmospheric 
hierarchy of splittings, Am^^^ ^ Am^^. 

In a recent paper [0, we performed a similar analysis in the SM framework (in 

which the only particles added to the SM are three right-handed neutrinos), also with 

positive results. There are similarities and differences between the SUSY and the SM 

cases. First, SUSY introduces additional unknowns in the scenario, particularly the 

supersymmetric mass spectrum. In the analysis, the only role of this spectrum is to 

give the threshold scale(s) below which the effective theory is just the SM. In this sense, 

the combination of (negative) experimental data and naturalness requires -^susy ^ 1 

^An alternative possibility to generate small non-zero neutrino masses in the MSSM involves im- 
parity breaking . See Ref . |^ for a discussion on the interplay between the two possibilities. 
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TeV. Of course, there may be appreciable differences between tlie masses of various 
supersymmetric particles (squarks, gluino, charginos, etc.). Still, the variation in the 
results of the RG analysis are not important. Hence, we will take a unique thresh- 
old at MsusY = 1 TeV throughout the paper. Let us mention here that, apart from 
three right-handed neutrinos, we will assume a minimal spectrum of particles; in other 
words we will work within the minimal supersymmetric standard model (MSSM). Sec- 
ond (and more important), in the supersymmetric regime, the charged- lepton Yukawa 
couplings are multiplied by a factor 1/ cos/3 with respect to their SM value. These 
couplings (together with the neutrino Yukawa couplings) play a major role in the ra- 
diative modification of the form of A4u- Thus, the results are going to present a strong 
dependence on tan j3 (we recall that tan /3 is defined as the ratio of the expectation 
values of the two supersymmetric Higgs doublets, tan/5 = {H2)/{H^))- Finally, the 
renormalization group equations (RGEs) themselves are different in the SUSY and in 
the SM cases. The difference is not just quantitative (i.e. differences in the size of the 
various coefficients), but also qualitative. In particular, the modification of the M.i, 
texture due to the contribution from the charged-lepton Yukawa couplings has oppo- 
site signs in the two cases (while the contribution from the neutrino Yukawa couplings 
themselves remains with the same sign). 

Let us briefly review the current relevant experimental constraints on neutrino 
masses and mixing angles (a more detailed account is given in ref. 0). Observations 
of atmospheric neutrinos are well described by — oscillations driven by a mass 
splitting and a mixing angle in the range [|1^ 

5 X 10"^ eV^ < ^m\^ < 10^^ eV^ , 

sin^ 29at > 0.82 . (1) 

Concerning the solar neutrino problem, as has been shown in ref. the small angle 
MSW solution is unplausible in a scenario of nearly degenerate neutrinos, so we are 
left with the large angle MSW (LAMSW) and the vacuum oscillation (VO) solutions, 
which require mass splittings and mixing angles in the following ranges 

LAMSW solution: 

10"^ eV^ < Aml^i < 2 x 10^^ eV^ 

0.5 < sin^ 29,oi < 1- (2) 
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VO solution: 

5 X 10"^^ eV^ < Aml^i < 1.1 x 10"^° eV^ 

8111^29^01 > 0.67. (3) 

From the previous equations, it is apparent the hierarchy of mass sphttings between 
the different species of neutrinos, Am^^; ^ Am^^, which should be reproduced by any 
natural explanation of those splittings. Let us also remark that it is not clear at the 
moment the value of the upper bound on sm^ 29soi [see eq.(@)]. As we will see, an 
upper limit like sin^2^so; < 0.99 or even greater, may disallow the scenario examined 
in this paper. For the moment, we will not consider any upper bound on sin^ 29soi (see 
ref.[y] for a more detailed discussion). On the other hand, according to the most recent 
combined analysis of SK + CHOOZ data (last paper of ref. [1^) the third independent 
angle, say (the one mixing the electron with the most split mass eigenstate), is 
constrained to have low values, sin^ 2(f) < 0.36 (0.64) at 90% (99%) C.L. 

Other relevant experimental information concerns the non-observation of neutrino- 
less double /3-decay, which requires the ee element of the matrix to be bounded 

as ini 

Mee < B = 0.2 eV. (4) 
In addition. Tritium /?-decay experiments indicate m^^. < 2.5 eV for any mass eigenstate 



with a significant Ue component |T2|. Finally, concerning the cosmological relevance of 
neutrinos, we will take J2 = 6 eV as a typical possibility and we will explain how 
the results vary when this value is changed. 

Let us introduce now some notation. In the SUSY framework the effective mass term 
for the three light (left-handed) neutrinos in the flavour basis is given by a term in the 
superpotential 

Weff = \v^M.V . (5) 

The mass matrix, M.y^ is diagonalized in the usual way, i.e. M^, = V*DV\ where 
D = diag(mie*'^, m2e*'^', ms) and is a unitary 'CKM' matrix, relating flavour to mass 
eigenstates 

C2C3 C2S3 
-C1S3 - SiS2C^e'^ C1C3 - SiS2S^e'^ 
S1S3 - CiSaCge*"^ -S1C3 - CiS2S'ie'^ 



S2e^'^ \ 






S1C2 


^2 ■ 


(6) 


C1C2 ) 







Here Si and Cj denote sin6'j and cos6i, respectively. In the following we will label the 
mass eigenvectors Ui as m^^ < m^^ and | Am^2 1 < I ^^23 1 ! where Amfj = mj 
is thus the most split eigenvalue). In this notation, constraint (H) reads 



2 A2S\ 



''1^2 ^2'^3 



V3 



(7) 



As it has been put forward by Georgi and Glashow in ref. 0, a scenario of nearly 
degenerate neutrinos should be close to a bimaximal mixing, which constrains the 
texture of the mass matrix TVI,^ to be essentially P,B[ 
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where niu gives an overall mass scale. M.h can be diagonalized by a y matrix 
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-1 
1 



V 



(9) 



2 2 v^y 

leading to exactly degenerate neutrinos, D = rriy diag(— 1, 1, 1), and 62 = 0, sin^ 26^3 = 
sin2 2^i = 1. 

It is quite conceivable that M.^ = M.h could be generated at some high scale by 



interactions obeying appropriate continuous or discrete symmetries [|T^. However, in 
order to be realistic, M.1, at low energy should be slightly different from M.^ to account 
for the mass splittings given in eqs.(|I|-^. We will explore whether the appropriate 
splittings (and mixing angles) can be generated or not through radiative corrections; 
more precisely, through the running of the RGEs from the high scale down to low 
energy. As discussed in ref.0, the output of this analysis can be of three types: 

i) All the mass splittings and mixing angles obtained from the RG running are in 
agreement with all experimental limits and constraints. 

a) Some (or all) mass splittings are much larger than the acceptable ranges. 

in) Some (or all) mass splittings are smaller than the acceptable ranges, and the rest 
is within. 
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Case (i) is fine, while case (ii) is disastrous. Case (in) is not fine, but it admits 
the possibility that other (not specified) effects could be responsible for the splittings. 
Concerning the mixing angles, it has been stressed in refs. |^,^ that due to the two 
degenerate eigenvalues of Aib, H is not uniquely defined. Hence, once the ambiguity 
is removed thanks to the small splittings coming from the RG running, the mixing 
angles may be very different from the desired ones. If such cases correspond to the 
previous (in) possibility, they could still be rescued since the modifications on A^,^ (of 
non-specified origin) needed to reproduce the correct mass splittings will also change 
dramatically the mixing angles. 

In section 2, we examine the general case in which the neutrino masses arise from an 
effective operator, remnant from new physics entering at a scale A. In this framework, 
we assume a bimaximal-mixing mass structure at the scale A as an initial condition 
and do not consider possible perturbations of that initial condition coming from the 
new physics entering at A. If tan f3 is small the LAMSW scenario in this case is of the 
undecidable type [possibility (Hi) above], but for tan/5 above a certain value which we 
compute, it is excluded (the VO solution is excluded for any tan/5). 

In section 3 we consider in detail a particularly well motivated example for the new 
physics beyond the scale A introduced before: the see-saw scenario. We include here 
the high energy effects of the new degrees of freedom above the scale A (identified now 
with the mass scale of the right-handed neutrinos) . We find regions of parameter space 
where the neutrino spectrum and mixing angles fall naturally in the pattern required 
to explain solar (LAMSW solution) and atmospheric neutrino anomalies, which we find 
remarkable. We complement the numerical results, presented in section 4, with compact 
analytical formulas which give a good description of them, and allow to understand 
the pattern of mass splittings and mixing angles obtained. We also present plausible 
textures for the neutrino Yukawa couplings leading to a good fit of the oscillation data. 

Section 5, still in the see-saw framework, discusses several possible implications of 
the effect of large neutrino couplings on: triviality limits on the Majorana mass; the 
infrared fixed point value of the top Yukawa coupling, with consequences for the lower 
limit on tan/5 and the value of the Higgs mass; and gauge coupling and bottom-tau 
unification. Finally we draw some conclusions. 
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2 M-jy as an effective operator 



In this section we will simply assume that the effective mass matrix for the left-handed 
neutrinos, J^^, is generated at some high energy scale, A, by some unspecified mech- 
anism. Assuming that below A the effective theory is the MSSM with unbroken 
i?— parity, the lowest dimension operator in the superpotential producing a mass of 



this kind is JTl 



W, 



4 2 2? 



(10) 



where k is a matricial coupling and if^ is the neutral component of the Y = +1/2 
Higgs field (the one coupled to the m— quarks). Obviously, A4,y = ^k{H2)'^. Between A 



and MsusYi the effective coupling k runs with the scale with a RGE [|14 
,dK 



dt 



+ 



K 



where t = log /i, and g2, Qi, Yt, Ye are the SU (2) x f/(l)y gauge couplings, the (MSSM) 
top Yukawa coupling and the (MSSM) matrix of Yukawa couplings for the charged 
leptons respectively. Between Msusy and Mz n runs with the SM RGE 



levr^ — 
dt 



-?>gl + 2A + 6y/ + 2TrYt Y 



K 

2 



kYIY,, + (YlY, 



K 



(12) 



where A is the SM quartic Higgs coupling and Yj, Ye correspond to the SM Yukawa 
couplings [the matching at Msusy is Yt{MsusY) = sin PYt{MsusY), Ye{MsusY) = 
cos /3Y,{Msusy)]. 

In a scenario of almost degenerate neutrinos, the simplest assumption for the initial 
form of the matricial coupling, k(A) is just the bimaximal mixing texture of eq.(|^), 
and this was also the assumption made in refs.[H,p[. In consequence. 



My 



(13) 



where = {Hlf + {H^f = (175 GeV)^. The last terms in eqs.(|lT], |T2D, i.e. those 
depending on Y|Ye, are generation-dependent and will modify the Ai^ texture, thus 
generating mass splittings and changing the mixing angles. It is easy to see MM that. 



in first approximation, the splittings Am?- = m^ — mf have the form 



Ami2 



--Am^ 
2 
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— AmL 
3 



-mte > 0. 



(14) 



where, neglecting for a while all the charged lepton Yukawa couplings but Yr, and 
working in the approximation of constant RGE f3 functions, e is given by 



~, 2 



Yr 



327r2 



, ^ , 1 MsusY 
log — h log 



(15) 



C0S2 p MsuSY /iO 

Thus, the SUSY and the SM corrections have opposite signs. If {K/ MsusyY^ ^^^"^ ^ > 
Msusy/ f^o, as it is the usual case, e has negative sign and the most split eigenvalue 
is the smallest one [thus the convention of labels used in eq.(|14|)]. Notice also that 
the pure SM case is recovered setting Msusy = A. As in the SM case, the previous 
spectrum is not realistic, i.e. the splittings are barely able to reproduce simultaneously 
the Aml^i and Am^^ splittings given in eqs.(|l]-|). Actually, the size of the splittings 
is larger than in the pure SM case (due to the coefficients of the RGEs and, especially, 
to the dependence on cos/3). This makes the scenario potentially more difficult than 
the SM one (see the discussion in the Introduction). 

Fig.l shows the complete numerical evaluation of the RGEs for mj, = 2 eV and 



tan/5 = 2, which corroborates the structure of eqs.([T^, |T5]). The splittings are always 
much larger than those required for the VO solution to the solar neutrino problem, 
Am^Q^ ~ 10~^° eV^. Therefore, the effect of the RGEs for this scenario is disastrous 
in the sense discussed in the Introduction for the possibility (ii). In consequence, as 
for the pure SM case, the VO solution to the solar neutrino problem is excluded^. For 
the LAMSW solution to the solar neutrino problem, things are a bit different. The 
smallest mass splitting is always within (or not far larger than) the LAMSW range, 
which indicates that we are in the situation {iii) explained in the Introduction (not 
satisfactory, but it may be rescued by extra physics). When is varied the results 
change with the scaling law Amfj oc m^. 

The dependence of the splittings on tan/5 is quite strong. Fig.2 shows this de- 
pendence for A = 10^° GeV. Even for moderate values of tan/5 the smallest splitting 
is much larger than the LAMSW range, thus spoiling that solution. Therefore, for a 
given value of A and m,^, the viability of the supersymmetric scenario of nearly degen- 
erate neutrinos puts an upper bound on tan/5. A reasonable estimate for this upper 
bound is given in Fig.3, which shows the value of tan/5 (vs. A) for which the small- 



^ Notice from the figure and from eq.(p^) that there is a value of A (close to Msusy) for which 
the splittings vanish, but of course the fine-tuning required for A to be close enough to that point, so 
that the splittings are within the VO range, is enormous. 
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A(GeV) 

Figure 1: Dependence of neutrino mass splittings at low energy (Am?- in eV^) with the cut-off scale 
A(GeV). For this figure tan/3 = 2. 

est splitting generated is larger (and ten times larger), than the maximum acceptable 
LAMSW value. The upper curve is therefore a conservative bound. These bounds 
become weaker if the neutrino mass, rrii,, decreases, as it follows from the previously 
mentioned scaling law. 

Concerning the mixing angles, in a first approximation the 'CKM' matrix, V, is 
given by 



/ 1 


1 


± \ 


V2 






1 




1 


2 




2^3 


1 







V 2 





(16) 



which leads to mixing angles 

sin2 2^i = ;|, sin2 2^2 = ^, sin2 2^3 = ^. (17) 
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tan (3 

Figure 2: Dependence of neutrino mass splittings at low energy {Amfj in eV^) with tan/?, for 
A = 10i°GeV. 

These values are far away from the bimaximal mixing ones. In consequence, they are 
not acceptable, stressing the fact that the simplified scenario discussed in this section 
does not work. However, as explained in the Introduction, if extra effects are able to 
modify the form of A^^ in order to produce the correct splittings, this modification 
will also change drastically the mixing angles, hopefully in a positive direction. This 
leads us to the see-saw scenario, which is analyzed in the next section. 



3 A4j^ from the see-saw mechanism 

The simplest example of the kind of new physics appearing at a scale A which can 
generate an effective mass term for the low-energy neutrinos we observe is the so-called 
see-saw mechanism [|^. Its super symmetric version has superpotential 

W = Wmssm - ^4-^4 + '^R^uL ■ H2, (18) 
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Figure 3: Upper limit on tan /3, as a function of A, beyond which the smallest neutrino mass splitting 
is larger (lower curve), or ten times larger (upper curve), than the maximum acceptable LAMSW value. 

where Wmssm is the superpotential of the MSSM. The extra terms involve three ad- 
ditional neutrino chiral fields (one per generation; indices are suppressed) not charged 
under the SM group: Ua^R {a = e,fi,T). Y^, is the matrix of neutrino Yukawa cou- 
plings and H2 is the hypercharge +1/2 Higgs doublet. Now, the Dirac mass matrix is 
niD = Yj,f sin /3. Finally, 7V1 is a 3 x 3 Major ana mass matrix which does not break 
the SM gauge symmetry. It is natural to assume that the overall scale of Ai, which we 
will denote by M, is much larger than the electroweak scale or any soft mass. Below 
M the theory is governed by an effective superpotential 

1 



Weff = Wmssm + ^{Y,L ■ H^fM-^Y.L ■ H2), 



(19) 



obtained by integrating out the heavy neutrino fields in (pISl). From this effective 
superpotential, the Lagrangian contains a mass term for the left-handed neutrinos: 



6C = --u'^Mjyiy + h.c., 



(20) 
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with 



Mu = uiu^M'^mj, = YjM'^Y,{Hl)\ (21) 

suppressed with respect to the typical fermion masses by the inverse power of the large 
scale M. 

The approximate degeneracy of neutrino masses in this framework follows from the 
initial condition for the neutrino mass matrix at the Planck scale: M.y{Mp) = Ai^, 
which would lead to exactly degenerate neutrino masses. We do not address in this 
paper the possible origin of this very symmetric form. As explained in 0, to explore 
the simplest (and most natural) textures of Yj, and Ai leading to this initial condition 
it is enough to consider the case in which the Majorana mass matrix is 



M = M 10, (22) 




and all the structure is in the Yukawa matrix, which reads 

Y, = Y,W,^. (23) 
Here Y^y is the overall magnitude of Yj^ and i? is a combination of two 'boosts' 



cosh a sinh a 
B= \ 1 

sinh a cosh a 



( cosh 6 sinh 6 
sinh 6 cosh 6 I , (24) 



V 1 

with two free parameters a, h. Having all the structure in the Majorana matrix and 
Yj, oc I3 is equivalent to the case a = 6 = 0. We will present our results, for fixed 
values of M and tan/5, in the (a, &) plane demanding |a|, |6| < 1.5. If \a\ or |6| are 
larger than 1.5 the matrix elements of Yj, are fine-tuned at least in a 10% 0]. 

The neutrino masses, exactly degenerate at tree-level by assumption, will receive 
generation dependent radiative corrections which will lift that degeneracy, eventually 
reproducing the pattern of mass splittings necessary to interpret the experimental in- 
dications. The bulk of these radiative corrections is logarithmic and easy to compute 
by standard renormalization group techniques: one starts at Mp with (^) and (^) as 
boundary conditions and integrates down in energy the relevant RGEs in a supersym- 
metric theory which is the MSSM with three right-handed neutrino chiral fields and 
the superpotential (|T^. At the scale M the v^^ are decoupled and below this scale 
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the running is performed exactly as in the previous section, except that now we have 
precise boundary conditions for the couphngs and masses at M. 

From Mp to M the evolution of the relevant matrices is governed by the following 



renormalization group equations 
dY„ 1 _ 



dt 



167r2 



3^2 + ^^1 - T2 



(25) 



rfYe 
dt 



167r 



2 Ye 



?>gl + ^-gl - Ti) I3 - (Y.tY, + SY+Y, 



where 



Ti = Tr(3Yl,YD + YtYe), T2 = Tr(3Y{jYu + YtY,), 



and 



dM 



M{Y,Yjf + Y,Y,^M 



(26) 



(27) 



(2^ 



(not yet given in the literature). Here g2 and gi are the SU{2)l and U{1)y gauge 
coupling constants, and Yu,D,e are the Yukawa matrices for up quarks, down quarks 
and charged leptons. 

At M, decouple, and Yg must be diagonalized to redefine the flavour basis of 
leptons [note that the last term in ( p6|) produces non-diagonal contributions to Yg] 
affecting the form of the Y^ matrix. Then the effective mass matrix for the light 
neutrinos is Mu ^ Y / M'^Y ^,{H^f . 

From M to M^, the effective mass matrix A1,y is run down in energy exactly as 
described in section 2. 

The renormalization group equations are integrated with the following boundary 
conditions: M. and Yy are chosen at Mp so as to satisfy 



My{Mp)=M,, 



(29) 



with the overall magnitude of Yj, fixed, for a given value of the Majorana mass M, by 
the requirement ~ (9(eV). We will take niy = 2 eV as a guiding example. The 
boundary conditions for the other Yukawa couplings are also fixed at the low energy 
side to give the observed fermion masses. The free parameters are therefore M, tan f3, a 
and b. 
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3.1 Analytical integration of the RGEs 



It is simple and very illuminating to integrate analytically the renormalization group 
equations (pT| , ^ , |26| , |28D in the approximation of constant right hand side. In this ap- 
proximation (which works very well for our analysis), the effective neutrino mass matrix 
at low-energy is simply Aii, plus some small perturbation. The overall mass scale is 
fixed to be of order 2 eV, so we need only to pay attention to the non-universal terms 
in the RGEs, which will be responsible for the mass splittings. Neglecting the Yg, 
Yukawa couplings, we get the following analytical expressions (the labelling of mass 
eigenvalues below may not always correspond to the conventional order ml < m 



U2 ' 



|Am22| < lAm^gl): 

m„ 1 + 3e^ - clclty ± 



tr + {clcl 



C2a 



+ 



1/2 



(30) 

where Ca = cosh a, S2a = sinh2a, etc. These expressions are identical to the ones 
derived in for the SM case, except for the numerical values of e,- and e,,, which are 
now given by: 



1287r2 



-2 log 



Mr, 



MsuSY 

Y^ 



+ cos^ (3 log 



M, 



SUSY 



— log — -. 

7r2 ^ M 



(31) 



(32) 



The Yukawa couplings in these expressions should be chosen at an appropriate inter- 
mediate scale but the simple formulas ( ^p2D are good enough for our purpose. Note 
also that Y^ and Yj. have an implicit dependence on tan f3 and, furthermore Y^, depends 
strongly on M (due to the requirement rriu ~ 2 eV). 

There are important differences with respect to the SM case presented in ||^: 1) 
is insensitive to the Majorana threshold, its sign is negative (it was positive in the 
SM) and it grows in magnitude for increasing tan/?; 2) e^, is twice larger than in the 
SM and decreases slightly when tan/5 increases. 

In the case a = b = 0, the mass splittings are independent of e^, and similar to those 
found in section 2, that is, not satisfactory (remember that this case is equivalent 
to having all the structure in A4, while Y,^ is proportional to the identity and thus 
universal) . 
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For a given M, there is a critical value of tan/3 below (above) which le,-! is smaller 
(larger) than e,y. This critical value of tan (3 increases with M. When ^ |eT-|, we can 
further expand the neutrino masses in powers of er/e^ finding 



Here we clearly see that the small mass splitting (solar) is controlled by the small 
parameter e,-, proportional to the squared Yukawa couplings of the charged leptons, 
while (proportional to the square of the larger Yukawa coupling Y^) is responsible for 
the larger mass difference (atmospheric). Moreover, the two neutrino eigenstates closest 
in mass (7711,-^,171^2) have masses of opposite sign, which is exactly what is required to 
fulfill the neutrinoless double /3-decay condition, eq.(0). In the SM case, for values of 
M between 10^ and 10^^ GeV, and e^y have the right orders of magnitude to account 
for AttiI^i, Am^j 0. In the supersymmetric scenario it is still the case that e,^ gives the 
correct atmospheric mass splitting, but for m,y = 2 eV gives typically a Aml^i one 
order of magnitude too high (recall here that by lowering mj, this splitting becomes 
smaller following the approximate law Am^ oc ml). This means that in this case 
only in particular regions of the (a, b) plane, in which some mild cancellation is taking 
place, we will obtain the right mass splitting to solve the solar neutrino problem. This 
cancellation does occur in regions around the lines a = and 2-\/2sinh6 = — tanha for 
which, in the approximation (|33|), m^^ = m^^- those regions then, we have a natural 
explanation for the Aml^i mass splitting while Am^j requires only a mild fine-tuning 
of parameters. 

If ^ |eT-|, then a different expansion shows that there is no natural hierarchy 
of mass splittings, which tend to be of the same order. In these stronger 
fine-tuning is required to get small enough Am^^j and the regions in parameter space 
where this occurs shrink, again around 1 — C2a — 2-\/2s2aS6 = 0. However, Am^j is still 
naturally of the right order of magnitude. 

Turning to the mixing angles, it is easy to see that the eigenvectors of the perturbed 



m, 



1 - + 4 + 




(33) 
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Aiu matrix are of the form 

Vl = Vi, vi 



1 



1 



-I3V2 + aV^ 



3), 



(34) 



where Vi are the eigenstates corresponding to the bimaximal mixing matrix, Vf, [see 
eq.®] 



(35) 
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In the approximation le,-! <^ Ci/, the parameters are given by 

CaSb + C(e,/e^) , /3 = + 0(e,/e^) . 



(36) 



The VI vectors define the new 'CKM' matrix V from which the mixing angles are 
extracted. In the case of interest, the spectrum of neutrinos consists of two lightest 
states iviivx-i) very close in mass and a heavier one rriy^. The relative ordering in mass 
of vtiy^ 2 is not fixed and this affects the relative order of Vl and i'^ ^^e V matrix. 
This has an effect on the sign of cos 26*3, which is important for the MSW condition 
cos 2^3 > (written using the conventional order m^^ < m^^)- This will be satisfied as 
long as V\ corresponds to the lightest mass eigenvalue. In other words, this condition 
requires that the negative mass eigenvalue, see eqs. (PO], PB| ), corresponds to the lightest 
neutrino. On the other hand, the relative ordering of the two lightest neutrinos does 
not change the values of sin^ 26*4 . 

In this approximation, if just one of the two (a, 6) parameters is vanishing, then 
V = Vb, i.e. exactly the bimaximal mixing case. Also, whenever c^, Cf, are sizeable 
(i.e. away from a = b = 0), \a\ ^ \(3\, and thus we are close to the bimaximal case. 
Therefore, it is not surprising that in most of the parameter space this will be in fact 
the case. This is remarkable, because it gives a natural origin for the bimaximal mixing, 
which was not guaranteed a priori due to the ambiguity in the diagonalization of the 
initial Ai u{Mp) = Aib matrix, as was explained in the Introduction. 

The only free parameter in the V matrix is the ratio a/ (3 and one obtains the 
relations: 



sin^ 29i 



(2r - 1)^ 
(2r + 1)^ 



sin^ 202 = 1 



sin^ 2^3 



1 - 



1 



2r)2' 



(37) 
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Figure 4: Upper plot: Neutrino mixing angles, sin^ 29i and sin^ 262, vs. sin^ 263 in the near degen- 
erate scenario. Lower plot: Zoom of the 0.9 < sin^ 29s < 1 region. 

where r = (3'^. It is instructive to invert the last relation, find r in terms of sin^ 29^, 
and then substitute back in the two previous relations. 

Figures 4a,b present sin^ 2 as functions of sin^ 29^ obtained in this way. For 
clarity, figure 4b focuses on the region of sin^ 29^ close to 1, which is the interesting one. 
We can impose the experimental hmits on these angles directly in this figure. We see 
that the hmit sin^ 2^2 < 0.64 (dashed fine) translates into a lower hmit sin^ 29^ > 0.94 
(the region to the left of the dashed fine is forbidden). In a similar way, sin^ 29 1 > 0.82 
(dotted fines) requires either sin^ 29s < 0.09 (but this solution gives sin^ 2^2 > 0.64 and 
is not acceptable) or sin^ 29^ > 0.9978 (the region to the right of the dotted line in figure 
4b is allowed). Note that if we impose the stringent condition sin^ 26*3 < 0.99, derived 
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Figure 5: Left plot: contours of Am^j/^V^ in the (6, a) plane from less than 10~^ (black area), 
through 6 x 10~^ (lines) to more than 2 x lO^"' (grey). Right plot: same for Am23/eV^, from 5 x 10~^ 
(black) to 10"^ (grey). The Majorana mass is 10^° GeV and tan/3 = 2. 

from some fits to solar data as discussed in tlie Introduction, then no viable solution 
would exist. Without such condition we find that this framework can accomodate the 
angles required by the data with 



sin^ 2^1 > 0.82, sin^ 202 < 0.17, sin^ 29^ > 0.9978. 



(3^ 



Since, in first approximation, the perturbed Ai^ will have eigenvectors of the form 
(3^), this result and the previous discussion are also valid for the SM case and any 
model starting with = M.^. Actually, the results shown in eg. (p7|) and Fig. 4 may 
be considered as predictions of such a kind of scenarios. 

In the next section we show that there are regions in parameter space where both 
the mass splittings and the mixing angles have the right values to explain the solar and 
atmospheric neutrino data. 
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Figure 6: Left plot: Contours of sin^ 262 in the (6, a) plane. The grey area marks the sin^ 26*2 > 0.64 
region. The line singled-out corresponds to sin^ 292 — 0.36. Right plot: Contours of sin^ 29i in the 
(6, a) plane. In the grey area sin^ 26i is smaller than 0.82, and the line corresponds to sin^ 26*1 = 0.9. 
The Majorana mass is 10^° GeV and tan/3 — 2. 

3.2 Numerical Results 

Figures 5 to 10 present our results for the mass splittings and mixing angles at low 
energy, after numerical integration of the RGEs from Mp to Mz as described in the 
previous sections. We stress that our convention in all the figures is to assign the indices 
1 and 2 to the neutrino states closest in mass, with m^^ < m^^ and leave the index 3 
for the other neutrino state. In the interesting regions, where the parameters would 
fit well the solar and atmospheric data, the third neutrino is the heaviest. Hence, the 
neutrino spectrum has the same structure as in the SM case 

We start with the results for M = 10^^ GeV and tan {3 = 2 chosen as a typical 
example; the dependence of the results with M and tan/5 is discussed afterwards. 
Figure 5, left plot, shows contour lines of constant Am^g (the squared mass difference 
between the lightest neutrinos) in the plane (6, a). The black (grey) region is excluded 
because there Am^2 < 10^^ eV^ (Amf2 > 2 x 10~^eV^), which is too small (large) to 
account for the oscillations of solar neutrinos (LAMSW solution). The white area is 
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thus the allowed region. The line in it corresponds to Amf2 = 6 x 10^^ eV^. Note 
how the allowed area is not too large and extends around the parametric line 1 — C2a — 
2\/2s2aSb = 0, as explained in the previous section. 

Figure 5, right plot, gives contour lines of constant Am2^. The black (grey) region 
is excluded because there Am^^ < 5 x 10~^eV^ (^^^23 > 10~^eV^), which is too 
small (large) to account for the oscillations of atmospheric neutrinos. Again, the white 
area is allowed. The small black areas correspond in fact to the "undecidable" case 
discussed in the Introduction: they might be rescued by unspecified extra effects. We 
do not present a plot for Am^g because it can always be inferred from Am^g and Am^2- 
Moreover, in the interesting case, Am^2 ^ Am23, one has Am^g ~ Am23. 

The intersection of the white areas in both plots is non-zero and would give the 
allowed area concerning mass splittings. It is always the case that the area surrounding 
the origin is excluded. There, the mass differences are always of the same order, and 
follow the same pattern discussed in section 2 (Am23 = 2Am^2)- case we 

conclude that, away from the origin, there is a non-zero region of parameter space 
where Am\^ ^ Amj2 5 accordance with the values required to explain the solar and 
atmospheric neutrino anomalies simultaneously. 

For mixing angles, figure 6, left plot, gives contours of constant sin^ 26*2 (one of the 
mixing angles relevant for atmospheric neutrino oscillations). The grey (white) area 
has sin^26'2 larger (smaller) than 0.64 and is disfavored (favored) by the data (SK + 
CHOOZ) at 99% CL. according to the most recent analysis (last paper of ref. ||1CI|| ). 
The line singled-out corresponds to sin^26'2 = 0.36 (maximum allowed value at 90% 
CL. according to the same reference). Figure 6, right plot, shows contours of constant 
sin^2^i (the other mixing angle relevant for atmospheric neutrinos). The grey (white) 
area corresponds to sin^ 20i smaller (larger) than 0.82, and is thus disallowed (allowed). 
The additional line included has sin^26'i = 0.9. 

Finally, figure 7, left plot, presents contours of constant sin^ 29^ which is relevant for 
oscillations of solar neutrinos. The grey (white) region has sin^ 26^ larger (smaller) than 
0.99. If one is willing to interpret the existing data as implying an upper bound of 0.99 
on sin^ 2^3, then the grey region would be excluded. The plotted curves give sin^ 26'3 = 
0.95. Figure 7, right plot, shows the region of the parameter space accomplishing the 
resonance condition (cos 2^3 > 0), which is required for an efficient MSW solution of 
the solar anomaly (see however the first paper of ref. [jl^] for caveats on this issue). 



19 



b b 




Figure 7: Left plot: Same as figure 6 for sin^ 263. The grey area corresponds to values above 0.99. 
The curves give sin^ 263 — 0.95. Right plot: The grey area corresponds to cos 2^3 < 

The region of parameter space where all constraints on mixing angles and mass 
splittings are satisfied is given by the intersection of all white areas in figures 5, 6 and 
7 (right plot). If sin^2^3 < 0.99 is imposed, then that intersection region, including 
now figure 7 (left plot), is empty and no allowed region remains. It should be noticed 
that this fact does not come from an incompatibility between the previous constraint 
and the sin^ 26*3 > 0.99 obtained from neutrinoless double /5-decay limits, eq. (0), in 
the 62 = approximation. If this were the case, it could be easily solved by decreasing 
the overall size of the neutrino masses, m^, in eq.(^, and this is not the case. Indeed, 
eq.(^ is satisfied in nearly all the parameter space. Even where sin^2^3 < 0.99, this is 
still true thanks to the contribution of 62- What actually forbids the whole parameter 
space if sin^26'3 < 0.99 is imposed is the incompatibility between acceptable 9i, 6*2 
and 6*3 angles to fit simultaneously all the neutrino oscillation data, as can be seen 
from the figures and in agreement with the discussion of the previous subsection, see 
eq. (p8D . This fact remains when rrii, is decreased. In fact, the effect of decreasing rriu is 
essentially an amplification of the figures shown here, which comes from the fact that 
for a given Majorana mass, the neutrino Yukawa couplings become smaller (the effect 
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Figure 8: Region (two disconnected parts) in the (6, a) parameter space for M = 10^° GeV and 
tan /J = 2 where all mass splittings and mixing angles satisfy experimental constraints. (See text for 
qualifications). 

is similar to decreasing M, which is discussed below). 

If the sin^2^3 < 0.99 condition is relaxed (as discussed in the Introduction), then 
the allowed region is given by the two islands in figure 8, which is non-negligible. It is 
remarkable that the small regions where all the mass splittings are of the right size and 
the small regions where the mixing angles are the appropriate ones have a non-zero 
overlapping. 

The dependence of the results with the Major ana mass M is illustrated by figure 
9, where the allowed regions in the plane {b, a) are shown for two different Majorana 
masses, M = 10^ and 10^^ GcV, keeping tan/3 = 2. For comparison with the allowed 
region for M — 10^° GeV (given in figure 8) we see that for small M the allowed 
region grows and files away from the origin, until it leaves the naturalness region 
|a|,|6| < 1.5. For M — 10^ GeV no allowed region inside the natural range for (a, 6) 
remains. Conversely, increasing M reduces the allowed region, which gets closer to the 
origin (at M — 10^^ GeV the allowed region becomes extinct). 

Concerning the remaining parameter, tan figure 10 gives the allowed regions for 



21 



1.5^^^^^^^^^^^^^^^ 1.5 




Q.5\ ^ 0.5 

(C 

-0.5^ I. -0.5 




-1.5 ^ \ -1.5 














-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 



Figure 9: 

Same as figure 8 for tan (3 — 2 and different values of tlie Majorana mass. Left plot: M = 10^ GeV; 
Right plot: 10" GeV. 

M = 10^° GeV and four different values of tan/3: 1.75, 3, 4 and 6.5, as indicated. 
The minimum in tan f3 is dictated by the requirement of perturbativity of all couplings 
up to the Planck scale, banning the presence of a Landau pole below it, as discussed 
in section 5. The allowed region gets smaller and smaller when tan/? increases. As 
explained in the previous section, \er\ grows with tan/5 making harder and harder a 
cancellation that gives the correct Am^^i and mixing angles. Eventually, for tan/3 > 6.5 
the allowed region disappears (for M = 10^^ GeV that value is tan/3 > 3.5). 

Finally, let us stress that, if the sin^2^3 < 0.99 condition is imposed, the whole 
parameter space becomes disallowed for any value of M and tan/3. We also find that, 
whenever there is a hierarchy in the mass splittings, the two lightest eigenvalues have 
opposite signs. This is just what is needed to have a cancellation occurring in the 
neutrinoless double /?-decay constraint (^. This constraint is satisfied in almost the 
whole parameter space for any M. 
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Figure 10: 

Same as figure 8 for M = 10^° GeV and different values of tan/3. Upper left: tan/3 = 1.75. Upper 
right: 3. Lower left: 4. Lower right: 6.5. 
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4 Examples of acceptable ansatze 



It is possible to find examples of matrices Y^{Mp) which both fall inside the allowed 
areas presented in the previous section and are natural (perhaps pointing to a possible 
underlying symmetry). 

An example, already presented in 0, which still survives in the supersymmetric 
case (for M ~ 10^° GeV and tan/5 ~ 2) is 



/ 



1 



Y,(Mp) = n 



2^2 
1 



V 



2^2 




1 



1 

1 1 

7^ 7t/ 



(39) 



It corresponds to a = and b = sinh ^(3/4) ~ 0.69. The mass splittings are 

Am?2 - 1 X 10"^ eV^ Am^g ~ 2 x 10"^ eV^ Am^g ^ 2 x 10"^ eV^ 
and the mixing angles 

sin^ 202 = 0.082, sin^ 29i = 0.9163, sin^ 29-^ = 0.99954, 



(40) 



(41) 



with cos 26*3 at the border of the resonance condition for the MSW mechanism. 

Another examples of working ansatze can be obtained. For instance, the following 
ansatz (corresponding to a = — cosh^"'^(v^/2) ~ —0.48, b = log(vT0/2) ~ 1.15 ) 

4 y/2 

1 Vb ^/5 



Y,(Mp) = n 



2^/5 

2 



(42) 



works correctly for M ~ 10^ GeV and tan /5 ~ 2, giving 

Am?2 ~ 1 X 10"^ eV^ Am^g ^ 9 x 10"^eV^ Aml^ - 
and the mixing angles 

sin^ 202 = 0.02, sin^ 201 = 0.979, sin^ 20^ = 0.99997, 



X 10"^ eV' 



(43) 



(44) 



with cos2^^3 > 0. It could be interesting to explore possible symmetries that may be 
responsible for the form of these ansatze and to analyze their implications for future 



long-baseline experiments |T8 . 
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5 Other relevant implications of neutrino-induced 
radiative corrections 



The previous sections focussed on the possibihty of reproducing the "observed" neu- 
trino mass sphttings from supersymmetric radiative corrections. In particular we ana- 
lyzed the case of nearly degenerate neutrinos, which is closely related to the bimaximal 
mixing scenario. 

In this section, still working in a see-saw framework, we take a different point of 
view and analyze the physical impact of neutrino-induced radiative corrections. The 
results of this section are quite generic. They are not associated to the bimaximal 
mixing scenario, not even to a scenario of degenerate neutrinos, but we will take this 
case as a representative example to illustrate the phenomena. Some of the following 
effects have already been mentioned along the paper, but here they are studied in 
greater detail. 

The first topic concerns the appearance of Landau poles associated to the neutrinos 
and the corresponding implications. It is easy to check from eq.(^) that the neutrino 
Yukawa couplings, Y,^, have a supersymmetric RGE quite similar to the top Yukawa 
coupling. It is therefore not surprising that they can develop Landau poles at high 
energy in a similar way. Obviously, the larger the right-handed neutrino Majorana 
mass, M, and the larger the low energy neutrino masses, m^, the larger the neutrino 
Yukawa couplings, thus lowering the scale at which the Landau pole appears. Con- 
sequently, if, in order to preserve perturbation theory (and the nice supersymmetric 
gauge coupling unification), we demand that the Landau pole does not appear below 
Mp (or the preferred high-energy scale), this puts upper bounds on M and m,^ (see 
ref.[0 for the analogue in the SM framework). 

To analyze this effect, we will use the simplest textures of Y^ and A^, leading to 
a bimaximal mass matrix, namely Y,^ = Y^l^ and M.{Mp) oc M.^ [M. has exactly 
the form given in eq.(|D but with overall scale M instead of m,^]. Recall that this is 
equivalent to the case a = 6 = in the analysis of textures performed in section 3. 
Obviously, for a, 6 7^ the neutrino Yukawa couplings are larger and the corresponding 
bounds stronger. Hence, we are analyzing here the most conservative case 

To extract the bounds, we set the Landau pole of at Mp (i.e. Yy{Mp) ^ 1) 
and evaluate the corresponding low energy value of niy, through the renormalization 
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M(GeV) 

Figure 11: Upper bound on the neutrino mass, ml^, vs. the Majorana mass M for two different 
values of tan/3. 

group equations of Y^, (between Mp and M) and k (below M) , for a certain value of the 
Majorana mass M. This "infrared fixed point" value, say m^^, represents an upper 
bound for the neutrino mass. The dependence of on M is illustrated in Fig. 11 
for two different values of tan/? (the value of M in the plot is to be understood as 
evaluated at the M— scale itself). Alternatively, for a particular value of the neutrino 
mass, nil,, we can extract from the figure the upper bound on the Majorana mass. 
We note that the bounds are quite strong. For quite moderate values of the neutrino 
masses, they conflict with the possibihty of a Majorana mass of 0{Mgut)- 

A different issue concerns the appearance of limits on tan/5. In section 2, it was shown 
that the scenario of nearly degenerate neutrinos is in conflict with large values of tan /3, 
the reason being that the radiative corrections to the mass splittings become much 
larger than the observed ones. The bounds were shown in Fig.3, for m,^ — 2 eV. For 
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Figure 12: Dependence of lowest value of tan/3 with the Majorana mass for = 2 eV. 



instance, supposing that the scale at which the effective mass operator is generated is 
A = 10^^ GeV, it can be seen from that figure that an upper limit tan (3 < 7 is obtained. 
This was fully confirmed by the see-saw analysis performed in section 3. As discussed 
in section 2, these upper bounds become weaker as the neutrino mass, m^, decreases. 

Now, from completely different reasons, neutrino-induced radiative corrections mod- 
ify the lower bound on tan /3. In the ordinary MSSM the lowest possible value for tan /3 
is derived from the value of the top Yukawa coupling in the infrared fixed point (IFP) 
limit, by the condition of a physical top mass in agreement with the experiment. In 
the MSSM extended with right-handed neutrinos, the RGE for becomes modified in 
the form 

dYu ^ Yu 

dt 167r2 

(we write the RGE for the full matrix of up-quark Yukawa couplings). If the neutrino 
Yukawa couplings, Y^, are sizeable (which occurs for large enough m,^ or M), their 



13 9 „ 9 16 9 ^ 

-9l + Sgl + Y9l-T2 



Is-SYIjYu-YIjYd 



(45) 
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contribution tends to lower the value of Yi further. In consequence, Ytimtop) decreases 
and the corresponding tan/5 increases. Fig. 12 shows the new lower bound on tan/5 
for the typical case my = 2 eV. Again, a diagonal structure for the neutrino Yukawa 
matrix, = Y^Is, has been assumed for simplicity. It is important to realize that we 
cannot raise the value of Y,^ arbitrarily, since then Y^ develops a Landau pole below 
Mp, as has been discussed in the previous subsection. This is the reason why the curve 
shown stops abruptly. As can be seen from the figure (evaluated at the 1-loop level), 
the change on tan/3 is very modest, but it may give rise to an increment on the Higgs 
mass of about 3 GeV, which should be taken into account for precision calculations. 

From a more qualitative point of view, we would like to mention two important im- 
plications of the presence of massive neutrinos for the supersymmetric perturbative 
unification. 

First, supersymmetric gauge coupling unification is considered as a brilliant success 
of the MSSM, since the two-loop running of the SU (3) x SU (2) x U{1) gauge couplings 
unifies with great precision at Mqut ~ 2 x 10^^ GeV. The unification, although re- 
markable, is not perfect. It is usual to invoke unknown (GUT or superstring) threshold 
effects in order to explain the discrepance. Since neutrino Yukawa couplings, Yj,, give 
a 2- loop contribution to the gi and g2 gauge couplings (see e.g. [|19|)7 for large enough 
Yjy (which means large my and/or M) this may be useful for a complete satisfactory 
unification. 

Second, the Y^ couplings have a 1-loop contribution to the charged lepton Yukawa 
couplings. Ye, see eq. (p6|) . Again, for large Y^, this produces significant variations for 
Ye at low energy, affecting the perturbative bottom-tau unification. This subject has 
been already addressed in the literature [^. Due to the form of eq. ([2B|) , working in 
the flavour basis for the charged leptons, the mass eigenvalues mg^ receive a correction 
from Yy that, at first order, goes like A^m^, ~ R ml_\YlYy]ii, where R is negative 
and flavour- independent. Thus, at this order, the neutrino contribution to the mi,/mr 
ratio goes always in the same sense for any Y^ texture. For the (problematic) first two 
families this means in particular that neutrinos are not useful to rescue the analogous 
second generation (m^/m^) ratio, since the correction goes in the wrong direction, but 
they might be useful for the m^/me ratio. In any case, non-trivial Y^ textures will 
have a significant impact on the supersymmetric mb/mr unification scenarios. 
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6 Conclusions 



We have studied the possibihty that neutrinos masses have 0{eV), in order to be 
cosmologically relevant. In that case they must be nearly degenerate, as required by 
the oscillation interpretation of atmospheric and solar data. An important question is 
whether radiative corrections have the right size to account for the small mass splittings 
required or are generically too large for the scenario to be considered natural. We 
addressed this problem in the context of the SM (plus three right-handed neutrinos) in 
a previous publication |^ and, in this paper, we have extended this thorough analysis 
to the supersymmetric case. 

The size of the mass splittings that we find is always much larger than required 
by the vacuum oscillation solution to the solar neutrino problem, solution which is 
therefore excluded in this scenario. 

When the origin of the non-zero neutrino masses is the see-saw mechanism (we 
concentrate our study in this appealing case) we find non-negligible regions in param- 
eter space where the mass splittings are consistent with the large angle MSW solution, 
providing a natural origin for the Am^^j <^ Am^^^ hierarchy. These regions correspond 
to Majorana masses around 10^'' GeV and small tan (3 ~ 2. Concerning the mixing 
angles, they are remarkably stable and close to the bimaximal mixing form (something 
that is not guaranteed a priory, due to an ambiguity in the diagonalization of the initial 
matrix) . 

We have understood analytically the origin of these remarkable features, giving 
explicit expressions for the mass splittings and the mixing angles. In particular, we 
give simple analytical relations between the mixing angles, which are also valid for any 
model starting with a bimaximal mixing form for the neutrinos. In addition, we have 
presented particularly simple see-saw ansatze consistent with all atmospheric and solar 
neutrino observations. 

Let us remark that the viability of the scenario is very sensitive to a possible up- 
per bound on sin^ 29^ (the angle responsible for the solar neutrino oscillations) . An 
upper bound such as sin^ 203 < 0.99 would disallow completely the scenario of nearly 
degenerate neutrinos due to the incompatibility between acceptable mixing angles to 
fit simultaneously all the neutrino oscillation data. This incompatibility is also easily 
understood from the analytical relations between the mixing angles, and thus is applies 
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also for the SM case and more generic models. 

Finally we have described in some detail several implications of the existence of 
(possibly large) neutrino Yukawa couplings. This includes the effects on: the trivial- 
ity limits on the see-saw Majorana mass, the infrared fixed-point of the top yukawa 
coupling, and the gauge and bottom-tau unification. 
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